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Abstract 

We study the geometry of billiard orbits on rectangular billiards. A truncated 
billiard orbit induces a partition of the rectangle into polygons. We prove that 
thirteen is a sharp upper bound for the number of different areas of these poly- 
gons. 
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1. Introduction 

Let a billiard ball be shot from a corner of a rectangular billiard. Consider 
the ball as a point, and truncate the orbit somewhere at the boundary. The 
truncated orbit of the ball generates a partition of the rectangular billiard into 
polygons, similar to Figure Q] Many of these triangles and quadrangles seem 
to have the same shape and size. In this paper we will show that (for a fixed 
shooting angle and stopping point) the number of different areas is at most 
thirteen. This universal upper bound is the sharpest possible. We also consider 
rational shooting angles and irrational shooting angles for which the thirteen is 
never reached. 

2. Rotations 

The results in this paper are closely related to the Three Gap Theorem (see 
e -g- 0; 0) an d the Four Gap Theorem (see [l[). The statements of these two 
theorems are best illustrated by a picture; see Figure [2] 

The Three Gap Theorem is naturally associated to the concept of rotations. 
First we recall the theorem and then we discuss rotations on intervals. For 
x € K, let {x} — x — [x\ denote its fractional part. 

Theorem 1. (The Three Gap Theorem) Let n €E N and a £ (0, 1). The num- 
bers 

0,W,{2a},{3a},...,{na} (1) 
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Figure 1: Truncated orbit of a billiard ball. The arrows indicate start and end of the orbit. 
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Figure 2: Left figure, the Three Gap Theorem: Cutting a pie n times where each next cut is 
obtained by shifting the previous one over a fixed angle a gives at most three different sizes 
of pieces of the pie. Right figure, the Four Gap Theorem: Now the first cut (at 0) works as a 
'reflecting boundary'. As soon as it is reached, we continue in the opposite direction. In this 
case we have after n cuts at most four different sizes. 
For this picture we used a = 0T405 * 2ir and n = 17. 

induce a partition of the interval [0, 1] in subintervals which can have at most 
three different lenghts. If there are three lengths, then the largest is the sum of 
the other two. 

Letting T a (x) — {x + a} for x € [0, 1], the numbers (JTJ) transform into 

o,r tt (o),r Q 2 (o),...,r; i (o). (2) 

If we consider i as a point on the circle of unit circumference, then T a {x) is 
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obtained by rotating x over a distance a. This gives a more dynamical view 
of the partition of [0,1]: the partition is induced by a truncated orbit of the 
rotation map T a . These observations lead to the following generalization of the 
Three Gap Theorem: 

Property 1. Let rii, n% € N, a € (0, 1) and a, b € R. The n\ + + 1 numbers 

aT- ni (0) + b, . . . , aT-\0) + b,b,aT a (0) + b,..., aT™ (0) + b (3) 

induce a partition of [b, b + a] in subintervals having at most three different 
lenghts. 

This can easily be obtained by taking n = n\ + in @, rotating over an 
appropriate angle and applying the linear map a ■ +b to the orbit. Actually, 
a special case of this property already appeared as a theorem in [l|. However, 
there a complicated proof was given to obtain this result. Vilmos Komornik 
came up with the idea to place the numbers on the circle, thus obtaining a 
much simplified and more natural argument [2j]. In the sequel we will refer to 
^ as a rotation orbit on [b, b + a] . 

There is a slightly stronger property we will need in Remark [Tj (see e.g. [l[ and 
1): 

Property 2. Take a truncated orbit of a rotation on an interval. Suppose the 
orbit consists of n numbers. Create another orbit from this by removing the last 
number. The two partitions induced by these orbits give two sets of lengths. The 
union of these two sets contains at most three different lengths. 

3. Billiards and the Four Gap Theorem 

The billiard in Figure [1] can be seen as a generalization to two dimensions 
of the pie-cutting process of the Four Gap Theorem, as illustrated in Figure [5J 
This statement deserves some explanation. Figure |3l a picture in some sense 
equivalent to the right panel of Figure [2l gives a description of the Four Gap 
Theorem in terms of a ball bouncing on the unit interval. 
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Figure 3: A ball bouncing on an interval between two walls. The Four Gap Theorem makes 
a statement about the subset of the interval consisting of the landing points of the ball. We 
used 0.1405 times the length of the interval as bouncing distance. 

This figure shows the movement of a ball bouncing between two walls, where 
we assume that the ball is a point and that there is no loss of energy. The landing 
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points of the ball build a sequence in the interval. The first n numbers in this 
sequence (0 included) define a splitting of the interval in n subintervals. The 
main statement of the Four Gap Theorem is that these subintervals can have 
at most four different lengths. In Figure [T] we now have a subset of a square, 
consisting of those points where the billiard ball appears. This observation gives 
already some reason to consider the billiard as a 2-dimensional generalization 
of the pie of the Four Gap Theorem. However, we can also argue this point of 
view in a more mathematical way. 

Let ||x|| denote the distance from x to the nearest integer. For a € M\ Q, let 

S a :=(\\k a \\)f =0 . 

Obviously this is a sequence in [0, i]. Moreover, it is exactly the sequence of 
landing points of a ball bouncing between and ^ with horizontal bouncing 
distance a. The sequence S a is obtained by 'folding' the sequence of integer 
multiples of a into the interval [0, h]. What we mean by this folding is illustrated 
in Figure SI where we plot the function fx : [0, oo) —> [0, |] 

fi{x) := \\x\\, 

and illustrate how [0, oo) is mapped to [0, ^] by f\. 




Figure 4: Plot of the folding map fi(x) = 



Now we concentrate on the billiard: the orbit of the billiard ball is obtained 
by 'folding' a halfline into a rectangle. Since the shooting angle is arbitrary 
between and 7r/2, we may without loss of generality assume that instead of 
a rectangle the billiard is a square and equal to [0, ^] 2 . The 'folding' map 
corresponding to this billiard is given by a two-variable function f2 : [0, oo) 2 — > 
[0, 1] 2 : 

/2(l,tf) = (|M|,||i/||). 

As we see, f2(x,y) = (fi(x), fi(y)), which is why the billiard can be viewed as 
being a generalization of the setting of the Four Gap Theorem to two dimensions. 
The folding map f% applied to a line creates a billiard orbit. Let a > 0, then 

B%, M] :={{\\x%\\ax\\):xe[Q,M]} 

describes a truncated billiard orbit that has initial slope a (the slope alternates 
between a and —a). Let A" M j and M j denote the number of different areas 
respectively different shapes in the partition of [0, i] 2 induced by M y Two 
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shapes are different if one can not be obtained from the other by translating, 
rotating and reflecting. For orbits truncated in a boundary point, we will prove 
the following theorem: 

Theorem 2. Let a > and choose M > such that (||M||, ||aM||) € [0, \] 2 \ 
(0, \) 2 ■ Then the billiard orbit B^ M j induces a partition of [0, in polygons 
for which 

Af 0>M] < 13 and 5[o,M] ^ 16 - 
These upper bounds are the best possible. 

In this theorem the billiard is square, but the result for rectangular billiards 
easily follows since the square can be scaled to any rectangle without changing 
the ratios between the shapes. From now on, we will assume that M satisfies 
the condition in the theorem. 

4. Orbit construction 

We already have an explicit expression for the billiard orbit B^ M j , but we 
will need a more tractable description. Therefore, in this section we present a 
rough intuitive outline of the way one can think of the geometry and construction 
of the billiard. The corresponding lemmata and their proofs are given in Section 
[H Consider the unit square and draw a line starting from the lower left corner 
with slope a. The boundaries are now considered to be connected as in a torus, 
so when we reach it, the line continues at the opposite boundary. Equivalently, 
if one of the coordinates is about to exceed 1, we subtract 1. But this is exactly 
taking fractional parts in both coordinates. Therefore, after we have traversed 
the unit square N times, we have a set which can be expressed as 

{{{x},{ax}) : < x < M} , 



for some Mel. A plot of such a set is shown in the left panel of Figure [5j 




Figure 5: Construction of a billiard orbit in three steps. Here N = 4 and a = v3 — 1. 

Now do the same starting from the other corners, traversing the square N 
times with a line either with slope a or —a. Explicit expressions for these four 
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sets (one for each corner) are given in Lemma [5J For an illustration, see the 
middle plot in Figure [5j 

The key observation now is that intersection of all AN lines with [0, i] 2 gives 
exactly a truncated billiard orbit with slope a, as is proved in Lemma 01 This 
fact is illustrated in the right plot in Figure [5j Obviously not all 47V lines actu- 
ally contribute to the billiard orbit. However, there is a good reason to consider 
them all: the intercepts of the 2JV lines with positive slope form a truncated 
orbit of a rotation on the interval [— a, 1], see Lemma El For the lines with 
negative slope a similar result holds. Having collected these insights, a simple 
counting argument suffices to obtain the upper bounds claimed in Theorem [2 

Proof of Theorem [2] Lemma [4] writes the billiard orbit as an intersection 
of the square [0, ^] 2 with a set of lines. Let us concentrate on the lines with 
positive slope. By Lemma |3] the intercepts of these lines form a rotation orbit 
on the interval [—a, 1]. So by Property [T] they induce a partition of this inter- 
val in subintervals of at most three different lengths. Denote the set of these 
lengths by D := {di, . . . ,d n }, where n < 3. For the lines with negative slope, 
the intercepts are the numbers 1 — yk, —N < k < N. They induce a partition of 
[0, 1 + a] in subintervals having lengths in the same set D. It now follows that 
vertical distances between adjacent parallel lines are in the set D. 




Figure 6: Local situation at the boundary where the orbit ends. Polygons of type 2 are 
triangular if the endpoint of the orbit is not one of the corners of the polygon (as is the case 
with Pi). There are only two shapes for which the endpoint of the orbit is one of the corners. 
One of them is still triangular (in this example P3), the other is irregular (P2). 

We will distinguish between three types of polygons: those that have no side 
which is part of the boundary of [0, ^] 2 (type 1), those that have exactly one 
such a side (type 2) and those that have two or more (type 3). 
The polygons of type 1 must be parallelograms. The area of such a parallelogram 
is given by didj/2a for some di, dj € D, and consequently they can have at most 
six different areas. 

A polygon of type 2 that is triangular must be half of a rhombus of which the 
vertical diagonal has length d £ D, and therefore its area is d 2 /4a. There is at 
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most one non-triangular type 2 polygon, as is explained in Figure[Hl So polygons 
of type 2 can have at most four different areas. 

Polygons of type 3 must be in one of the corners of [0, i] 2 , but not in (0, 0) since 
the orbit starts there. So this gives at most three more areas. 
Putting everything together, it turns out that the number of different areas is 
bounded by thirteen. 

For the number of shapes a similar counting argument holds. The number of 
parallelogram shapes is again six, since reflections do not count. The triangles 
that are half of a rhombus can have at most six different shapes, since there are 
three types of rhombi which can be cut either horizontally or vertically. The 
rest of the argument doesn't change, so there are at most three more different 
shapes than different areas, which establishes the upper bound of at most sixteen 
different shapes. 

The sharpness of these bounds follows from Example [T] in section [7] □ 

Remark 1. As the careful reader may have noted, the construction of the bil- 
liard orbit always gives a truncation on the left boundary or on the lower bound- 
ary of the square. So strictly speaking, Theorem[2]is not proved in full generality 
yet. Suppose we have an orbit truncated at the upper or right boundary. By 
removing the last linear part or adding the next linear part, we can transform 
this orbit into an orbit truncated at the left or lower boundary. This means 
that in the proof above, the rotation orbit on the interval [—a, 1] contains one 
element more or one less than the rotation orbit on [0, f + a]. Now Property [2] 
tells us that vertical distances between adjacent parallel lines can still have at 
most three different values, completing the proof. 

5. Rational angles and a golden exception 

Theorem [5] gives an upper bound for the number of different areas of shapes 
on the billiard table. Some natural questions remain. For example, what hap- 
pens if a is rational? Can we prove sharper upper bounds under suitable con- 
ditions? In this section we explore these properties. 

Obviously, taking a rational gives a special case. The first thing to note is that 
the orbit will be periodic: if a — p/q, then for x G R 

{\\x + q\\,\\a{x + q)\\) = {\\x\\, \\ax\\). 

A bit less trivial is the following result. 

Proposition 1. The best upper bound for A" Q M j with a G Q is 13, but for all 

a G Q there is an Mq such that 1 < .4^ M j < 3 for M > Mq ■ These bounds are 
sharp. 

Proof Note that the areas of the polygons continuously depend on a. So if we 
have an a and M such that A^ Q M j = 13, then we can find e > such that the 
upper bound of thirteen is reached for all a G (a — e, a + e). Since this interval 
contains rationals, we see that rationality is not sufficient for a sharper upper 
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bound. 

Since the orbit is periodic, the partition doesn't change anymore if M is large 
enough. Taking a = 1 shows that 1 is a sharp lower bound for the limiting 
number of shapes. For the upper bound, suppose that a = p/q. By Lemma [3] 
the intercepts satisfy 

Vp+ q = ( 1 + - )\ — — — : — \ - - = — = yo- 
V qj I l+p/q J q q 

It follows that the numbers yk form a periodic rotation orbit on [—a, 1] and 
therefore the set D as defined in the proof of Theorem [2] contains only one 
length if M is large enough. If p and q are relative prime, then this length 
is 1/q. Now a type 1 polygon is a rhombus with area l/2pq. Since there is 
no endpoint of the orbit anymore, a type 2 polygon is half of such a rhombus. 
Polygons in the corners are also triangular, because the orbit touches all sides 
of the square before becoming periodic. These triangles are quarters of the 
rhombus, thus having area l/8pq. This makes at most three different areas in 
total. To see that this upper bound is sharp, see Figure [7J □ 




Figure 7: The periodic orbit for a = 3/5. There are three different areas: the rhombi have 
area 1/(2 ■ 3 ■ 5) = 1/30. The triangles have area 1/60 or 1/120. 

Surprisingly, there exist irrational a for which the upper bound of thirteen 
different areas is never reached: 

Proposition 2. Let 4> = ("v5— l)/2 denote the small golden mean. If a = 
for some neN, then A? M i < 12. 

Proof Consider the numbers yk that form a rotation orbit on [—a, 1). The 
partition of [—a, 1] induced by this orbit gives subintervals with lengths in a set 
D. This set D changes if we extend the orbit (i.e. we increase M): some lengths 
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will disappear and new lengths will be created. In [1| and |3| it was shown that 
the largest length is always the first to disappear. A new length only pops up 
if there are only two lengths in D, and the new length is the difference of these 
two existing lengths. Together with the fact that 1 — <fi = <fi 2 , this is the basis 
of our argument. 

Let a — l/(n + 4>). From the way points are added to the rotation orbit it 
is clear that we can choose M such that [—a, 1] will be partitioned in n + 1 
intervals of length a and an interval of length 1 + a — (n + l)a = 4>a. This 
gives D = {a, (f>a}. Extending the orbit with one more point transforms D into 
{a, (j)a, </> 2 a} and this is the first time that D contains three lengths. Increasing 
M further, D will change into {<fia, 4> 2 a} and then into {4>a, 4> 2 a, </> 3 a} . An 
inductive argument suffices to show that the ratios between the lengths in D 
are preserved. 

Recall that the areas of the parallelograms are determined by a product of two 
lengths in D. By the above reasoning, if D = {di, c?2, ^3}, then did^ = d 2 ,, 
which implies that the parallelograms can have at most five different areas. 
Consequently A" M j < 12. □ 

6. Lemmata and their proofs 

Let a > be an irrational number and consider the halfhne l(x) — ax, x > 
0. Let Si — [0, l) 2 and define S%, S3, 1S4, . . . to be the squares of the form 
[k, k + 1) x [m, m + 1), with k and m integers, that are consecutively traversed 
by the halfline, see Figure [8] Choosing an index N, there exists M € R such 
that 

N 

[J S k n{(x,ax) : x > 0} = {(x,ax) : < x < M} . 

k=l 

Taking fractional parts in both coordinates can be seen as mapping each of 
the squares Sf. to [0, l) 2 . Therefore, doing this for the above set gives 

N 

{({x} , {ax}) : < x < M} = [0, l) 2 n [j {{x, ax + y k ) : x e K} (4) 

fc=i 

for numbers yk defined by the recursion 
Vi = 0, 

[ Vk + a if y k < I -a, . , 

[ yk-i if yk > 1 - a. 

We will denote the set in Q by A ++ . The ++ superscript reflects the fact 
that we started with a halfline in the first quadrant, so both coordinates are 
positive. Doing similar operations to halflines in the second, third and fourth 
quadrant, we can define sets A A and A^ respectively as follows: 

A-+ = {(I - {x} , {ax}) :0<x< M} , 
A— = {(1 - {x} , 1 - {ax}) : < x < M} , 
A+- = {({x},l- {ax}) :0< x < M} . 
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1 2 3 4 5 



Figure 8: Construction of the squares S fe . Here a = v3— 1, N = 8 and M = 5. The numbers 
yk are approximately given by yg = 0, j/i ~ 0.732, ?/2 ~ —0.268, ys ~ 0.464,. . . Compare with 
Figure [5] left plot. 



Taking the union of these four sets and intersecting with [0, ^] gives us a 
billiard orbit, as is proved in the lemma below. 



Lemma 1. The billiard orbit B? n , satisfies 



[0,M) 

1 2 



B fo,M) = U ^ n t°.^ 



2 J 

u,ue{+,-} 

Proof Observe that 

(\\x\\, ||ckx||) = ^min|{a;} , 1 — min|{aa;} , 1 — {aijj 

= [0,^] 2 n |J |J (a, 6), 

a£{{i),l-{i}} b£{{ax}, l-{ax}} 

and now take the union over all x £ [0, M). □ 

In the next lemma expressions similar to ((4]) are derived for A ^ , A and 
A+~. 

Lemma 2. Let y^k = 1 — ol — Vk for k = 1, 2, . . . , N. Then 

N 

A-+ = (0,1] x [0,1) n (J {(x,-ax + l-y„ k ) : x e K} , 

fc=i 
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N 

A— = (0,l} 2 n{J{{x,ax + y- k ):xeR} 7 

k=l 

N 

A+- = [0,1) x (0,1] n (J {(x,-ax + l-y k ) :xeR}, 

fc=i 

Proof Define the functions f,g, h : R 2 -» E 2 by f((x,y)) = (1 — x,y), 
g{{x, y)) = (1 — x, 1 — y) and y)) = (x, 1 — y). Applying these functions to 
the left hand side of Q, we get = ,g(A++) = A and h{A++) = 

A^ . On the other hand, 

f({(x,ax + y k ) :x eR}) = {(1 - x, ax + y k ) : x e R} 

= {(x,a(l -x) +y k ) : x e M} 
= {{x. -ax + 1 - y- k ) : x e K} , 

whence application of / to the right hand side of ((4]) leads to 

N 

/([0,1) 2 n [J {(x,ax + y k ) 
fc=i 

N 

= /([0,1) 2 ) n/(U {(x 7 ax + y k ):xeR}) 

k=l 
N 

= (0, 1] x [0, 1) n |J /({(x, ax + y k ) : x G R}) 

fc=i 

AT 

= (0, 1] x [0, 1) n [J {(x, -«x + 1 - y_ fe ) : x G M} , 
fc=i 

so for A *~ we established the equality claimed in the lemma. The other two 
equalities for A and A^ follow from a similar reasoning since 

g({{x,ax + y k ) : x € E}) = {(1 - x, 1 - ax - y k ) ■ x e R} 

= {(x, 1 - a(l - x) - y k ) : x € E} 
= {(x, ax + y~k) '■ x € M} , 

and 

ax + ?/fe) : x G E}) = {(x, I — ax — y k ) : x G E} . 

□ 

The numbers y k and satisfy a nice relation, as is shown in the following 
lemma. 

Lemma 3. Let yo = —a. Then the numbers y~N, • ■ • , Vn form a rotation orbit 
on the interval [—a, 1]. They are given by 

( ka 1 

y k = {l + a){-— \-a for -N<k<N. (6) 
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Proof The recursion ([5]) can be rewritten as 

Vk+i = (Vk + 2a mod(l + a)) - a, 

and therefore 

y k+ i + a y k + 2a ( y k + 2a 

mod 1 — 



1 + a 1 + a { I + a 

Letting y k = , k = —N,...,N and a — j^, this reduces to 

Vk+i = {yk + a} . 
Since y\ = 0, we have y\ = a, which leads to 

Vk = {ka} for k > 1. 
On the other hand, for k > 1, 

_ y-k + a _ 1 - a - y k + a _ l + a _ y k + a 
l + a l + a l + a l + a 

= 1 — 2/fc = 1 - {^a} = {-/ca} , 



since a is irrational. By definition we have yo = 0, and hence 

y k = {ka} for - N < k < N. 
Solving for y k gives the result. □ 

In Lemma [T] we already derived an expression for M ^ , but this is not so 
easy to analyze directly. In the next lemma we describe B^ M j as the union of 

two collections of lines intersected with [0, ^] 2 . All lines in the first collection 
have slope a and all lines in the second collection have slope —a. 

Lemma 4. Let lt(x) — ax + y k and l k (x) = —ax + 1 — y k . Then 

1 N 
£fo,M] = [o, n |J |J {(x,l%(x):xeK)} 

«£{+,-} k=-N 

Proof This lemma will be proved by taking closures in the equation in Lemma 

m 

°[<d,M) ~ -°[0,M]> 

since x h-> (||x||, ||ax||) is a continuous function from K to R 2 . On the other 
hand, 

[o,i] 2 n |J A- = [o,i] 2 n |J A™>, 

u,v£{+,-} u,v£{ + ,-} 
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and since A uv is a finite collection of lines intersected by a 'half open' unit square 
its closure is the same collection of lines but now intersected by the closed square 
[0,1] 2 . Therefore, 



JY 



A++ U A = [0, l] 2 n (J {{x,lt{x)) (7) 
k = -N 

Now note that since Iq(x) — ax — a we have 

[0,\} 2 n{(x,l+(x)):xeR}=<D. 



Intersecting both sides of ([7]) with [0, i] 2 gives 

N 

, 2 n 

2 J 2 



i i N 

[0,-] 2 nA++UA- = [0,-] 2 n |J {(x,l+(x)):xGR} (8) 

k=-N 

Analogously it follows that 

i i N 

[0,-] 2 nA+-UA-+ = [0,-] 2 n |J {(x,l^(x)):xER} (9) 

k=-N 

Combination of the last two equations gives the result. □ 



7. Sharpness of the bounds 

In this section we present an example in which the upper bounds of Theorem 
[2] are reached. This proves sharpness of the bounds. 

Example 1. Let a = and choose N = 11. The corresponding orbit is 
shown in Figure [9] Use Lemma [3] to find the numbers yk and let 

di = y-2 - Vi w 0.0965, d 2 = y n - y_ 2 » 0.0658, d 3 = y_ 5 - Vn « 0.0307 

denote the three different vertical distances between adjacent parallel lines. The 
areas of the shapes are of the following form: 

Shapes I, II, III, IV, V, VI : d i d j /2a, i<je {1,2,3} 

Shapes VII, VIII, IX : dj/Aa, ie {1,2,3} , , 

Shape X : d z d 1 /2a - df/ia [ } 

Shapes XI, XII, XIII : d 2 /8a, ie {1,2,3} 

Calculating these thirteen areas indeed gives thirteen different values, where 
a precision of two decimals suffices. The flakes VII, VIII and IX have the 
same areas as Vila, Villa and IX a respectively, so the maximal number of 
sixteen different shapes is also reached. We checked the calculations by using 
the outcomes to determine the area of [0, i] 2 . 
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Figure 9: Thirteen different areas, sixteen different shapes. 
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